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We study charge pumping when a combination of static potentials and potentials oscillating with a 
time period T is applied in a one-dimensional system of non-interacting electrons. We consider both 
an infinite system using the Dirac equation in the continuum approximation, and a periodic ring with 
a finite number of sites using the tight-binding model. The infinite system is taken to be coupled to 
reservoirs on the two sides which are at the same chemical potential and temperature. We consider a 
model in which oscillating potentials help the electrons to access a transmission resonance produced 
by the static potentials, and show that non-adiabatic pumping violates the simple sin0 rule which 
is obeyed by adiabatic two-site pumping. For the ring, we do not introduce any reservoirs, and we 
present a method for calculating the current averaged over an infinite time using the time evolution 
operator U(T) assuming a purely Hamiltonian evolution. We analytically show that the averaged 
current is zero if the Hamiltonian is real and time reversal invariant. Numerical studies indicate 
another interesting result, namely, that the integrated current is zero for any time-dependence of 
the potential if it is applied to only one site. Finally we study the effects of pumping at two sites on 
a ring at resonant and non-resonant frequencies, and show that the pumped current has different 
dependences on the pumping amplitude in the two cases. 

PACS numbers: 73.23.-b, 73.63.Nm, 72.10.Bg 



I. INTRODUCTION 

The idea that oscillating potentials applied to cer- 
tain points in a one-dimensional system can pump a 
net charge between two reservoirs at the same chemical 
potential has been studied extensively for many years, 
both theoretically [1-42] and experimentally [43M49j . For 
the case of non-interacting electrons, theoretical studies 
of this phenomenon have used adiabatic scattering the- 
ory |9rlr7|. Floquet scattering theory 21 23], the non- 
equilibrium Green function formalism [2414271 ] . and the 
equation of motion approach [3(| H3] • The case of inter- 
acting electrons has been studied using a renormalization 
group method for weak interactions [50j , and the method 
of bosonization for arbitrary interactions [51-61]. Other 
interesting studies of pumping include adiabatic quantum 
pumping in graphene where the electrons obey the Dirac 
equation [621 ] . and classical pumping on a finite ring by 
oscillating hopping rates at two sites [63l ]. 

With the exception of a few papers (22l. I24H271 Hfj j . the 
earlier studies of charge pumping have generally consid- 
ered systems in which oscillating potentials are applied to 
two or more sites. In such cases, it is known that if the os- 
cillation frequency w is small, the dc part of the pumped 
current is proportional to uj; the charge pumped per cy- 
cle (with time period 2tt/uj) therefore has a finite value in 
the adiabatic limit ui — > 0. However, it has been noted in 
Refs. [22| and I24H271 that an oscillating potential applied 
to a single site can also pump charge provided that the 
system has no left-right symmetry; this can happen if, for 
instance, appropriate static potentials are present. Most 
studies of charge pumping have also been limited to infi- 
nite systems in which the left and right sides of the system 
(called the leads or reservoirs) are associated with certain 



chemical potentials and temperatures. Charge pumping 
on a finite ring has been studied in a few papers for adi- 
abatic [lli - flfjj and nonadiabatic situations [5^, (2|| [26[ . 

In this paper, we will study charge pumping in both 
an infinite system as well as on a finite ring for non- 
interacting electrons. We will study the effects of oscil- 
lating potentials applied to either one site or more than 
one site. We will not assume the oscillation frequency 
to be small (i.e., the adiabatic limit), but will assume 
the oscillation amplitudes to be small. For the infinite 
system, we will assume the existence of reservoirs, but 
for the finite ring, we will assume that the system is not 
coupled to any reservoirs and that the time evolution is 
purely Hamiltonian. We will see that results obtained in 
the two cases differ in some interesting ways. We will re- 
strict our analysis to zero temperature and spinless elec- 
trons, since spin does not play an essential role in the 
absence of interactions between the electrons. 

The plan of the paper is as follows. In Sec. II, we 
will examine charge pumping in the infinite system. For 
convenience, this will be studied in the continuum limit 
using a linearized form of the energy-momentum relation; 
namely, we will use the massless Dirac equation with both 
right and left moving modes. We will assume the system 
to be coupled to reservoirs on the two sides, and the 
chemical potentials and temperature of these will be in- 
troduced using the formalism of Refs. and This 
implicitly assumes that an electron, after passing through 
the region with the static and oscillating potentials, equi- 
librates with whichever lead it enters; the precise mech- 
anism for energy or momentum relaxation in the reser- 
voirs will not be specified in our calculation. As a specific 
example, we will consider a model in which the static po- 
tentials have transmission resonances at certain energies, 
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and potentials oscillating at exactly the right frequency 
can then lead to enhanced charge pumping [38l . |3^ | ; in 
this model, we show how a competition between differ- 
ent processes can lead to either maxima or minima in the 
pumped current. We find that the pumped charge is of 
second order in the strengths of the oscillating potentials 
and can be much larger for two-site pumping compared 
to one-site pumping. We also find that the 'sin rule, 
which has been discussed earlier in the adiabatic limit 
uj — > for pumping by two oscillating potentials with a 
phase difference of <fi [H EH > breaks down if u is larger 
than the resonance width. 



In Sec. Ill, we will study charge pumping on a fi- 
nite ring using a tight-binding Hamiltonian to describe 
the electrons; no reservoirs will be introduced. We will 
present a formalism for calculating the long-time aver- 
aged current if there are potentials oscillating with a time 
period T; this is done by considering the eigenstates of 
the unitary time evolution operator U(T). Next, we will 
assume that the system begins in a specific initial state 
corresponding to a particular filling at zero temperature. 
We will then evolve the system using only the Hamilto- 
nian, without introducing any mechanisms for momen- 
tum relaxation and phase decoherence; this means that 
the system will never reach a steady state and that its 
properties depend on the initial state. We will show ana- 
lytically that if the Hamiltonian is real and time reversal 
invariant, then charge pumping cannot occur even if the 
system has no left-right symmetry. Numerically, we also 
find that if the Hamiltonian is real and if an oscillating 
potential is applied to only one site, then charge pumping 
does not occur even if the oscillating potential is not time 
reversal invariant and the system has no left-right sym- 
metry. It is therefore necessary to apply oscillating po- 
tentials to at least two sites in order to pump charge. We 
will also study what happens at both non-resonant and 
resonant frequencies; the latter means that the oscillation 
frequency is equal to the energy difference between a filled 
state and an empty state of the time-independent part of 
the Hamiltonian. (We would like to mention here that 
charge pumping on a ring at resonant and non-resonant 
frequencies has also been studied earlier in Ref. [23h . In 
the Appendix, we will present details of the calculation 
of the eigenstates of U(T) to first order in the oscillat- 
ing potentials, both for the non-resonant and resonant 
cases. We will show there that the pumped charge can 
receive contributions at either zero-th or first order in 
the oscillating potential in the resonant case but only at 
second order in the non-resonant case. In Sec. IV, we 
will summarize our results and emphasize the different 
assumptions that we have made in the models defined on 
the infinite line and on a ring. Finally, we will discuss 
possibilities for experimentally testing our results in two 
kinds of systems, namely, mesoscopic rings and aromatic 
molecules. 



II. CHARGE PUMPING ON AN INFINITE 
LINE 

In this section, we will study charge pumping on an in- 
finite line. Our model consists of a gapless system of non- 
interacting spinless electrons subject to some static and 
oscillating potentials. We will assume that the regions 
lying far to the left and far to the right of all the poten- 
tials have the same chemical potential given by the Fermi 
energy Ep (we will work at zero temperature). Assum- 
ing that the oscillating frequency is small compared to 
the bandwidth, we will linearize the energy-momentum 
dispersion around Ep. Let us denote the Fermi veloc- 
ity and Fermi wave number by vp — [dE j dp) e=e f and 
kp respectively. [The relation between hp and Ep is 
governed by the underlying microscopic model. For in- 
stance, in a tight-binding lattice model with the disper- 
sion E = — 27COsfc, where 7 is the nearest-neighbor 
hopping amplitude, we have Ep = — 2^/ cos kp where 
< kp < it]. We can then define the electron field 
operator 

*P(x) = 1> R (x) e ikpx + 4> L (x) e~ lkFX , (1) 

where %pL and ipR are the fermionic field operators for 
the left and right moving electrons. In terms of these 
fields, the Hamiltonian in the presence of several point- 
like potentials is given by H = Hq + V, where 

IT f J ■ I I t 9 ^ R , / t d ^L , 

H Q = J dx ivp (- 4> R -^r + tL-foh 

V = I dx 5{x-x p ) U p (t) ^{x)ijj{x). (2) 
J v 

In the absence of the potential V, the eigenstates of H are 
given by exp[i(±kx — Et)\, where E — vpk, and ±k refer 
to right and left moving modes respectively; note that we 
are measuring the energy E and the wave number k with 
respect to Ep and kp respectively, and we have set K = 1. 
Using Eq. (TTJ), we find that the equations of motion for 
ipn and ipL are given by [6l| 

. dip R . dtp R 



p 



= ]T 6(x-x p ) U p (t) (Vl + ^ fl e <a ***'). (3) 
p 

We can solve these equations by integrating over small 
regions from x v — e to x v + e to find the discontinuities 
in the fields at x — x p . 

Let us now consider two types of point-like potentials: 
a m S(x — x m ) which are time independent, and w n (t)d(x — 
x n ) which oscillate in time as w n (t) = b n cos(ujt + cf> n ). 
We will assume that the oscillation amplitudes b n are 
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small compared to vf and will generally calculate the 
pumped charge to the lowest non-zero order in the 6 n 's. 
Let us first assume that b n — for all n, and that we can 
completely solve the problem with the static potentials. 
This can be described in terms of a scattering matrix S 
as follows. If TZ denotes the region within which all the 
static potentials are present, then an electron incoming 
from a region far to the left of TZ (which we denote as x <C 
TZ) with unit amplitude and energy and wave number 
given by E and k = Ejvp respectively will have a wave 
function %j){x)e~ t given by 



= t R {E) e 



r L (E) e~ lkx for x « TZ, 
lkx for x^TZ, (4) 



where rp. and t R denote the reflection and transmission 
amplitudes. Similarly, an electron incoming from a region 
far to the right of TZ with unit amplitude and energy and 
wave number given by E and — k = —E/vp respectively 
has the wave function ip(x)e~ iEt , where 



—ikx 



ip(x) = e 

= t L (E) e 



- r R {E) e lkx for x > TZ, 

~ ikx for x<.TZ. (5) 



Assuming that all these reflection and transmission am- 
plitudes are known, let us now turn on the oscillating 
potentials, all of which we take to lie to the left of TZ, 
and proceed as follows. 

An electron incoming from the left of TZ with unit am- 
plitude and energy (wave number) Eq (kg = Eq/vf) will 
now have a wave function ^(x, t) of the form 



ip = e 



i(k x-E t) 



rL >3 e l{ ~ k]X ' E]t) for x < TZ, 



t RJ e i! - k i x - E ^ for x > TZ, 



(6) 



where Ej = Eq +joj, kj = Ej/vf, and tl,j and t R j will 
now generally be functions of both Eq and Ej . Similarly, 
an electron incoming from the right of TZ with unit am- 
plitude and energy (wave number) Eq (ko = Eq/vf) will 
have a wave function ip( x > t) of the form 

tJj = e i(-kox-E t) + r RJ e iQ>i*-Eit) for x » TZ, 

3 

= tL,ie i{ - kiX ~ Bit) for x < TZ, (7) 



where Ej = Eq + ju, kj — Ej/vF, and r R j and tL,j will 
be functions of both Eq and Ej . 

The sums over the side band index j in Eqs. ([6][7]) will 
go from — oo to oo for a Dirac electron. But, in practice, 
j gets cut off for two reasons. First, if the energy Ej goes 
above or below the bandwidth of the system, the corre- 
sponding wave function decays exponentially as \x\ — > oo; 
such a state does not carry any current. Secondly, if all 
the pumping amplitudes b n are small, one can show, us- 
ing Eqs. ([3]) recursively, that the leading order terms in 



r R /L,j and t R/LJ are given by |&„|I J I if j ^ 0; this goes 
to zero exponentially as \j\ — > oo. Hence, if w is much 
smaller than the bandwidth and Eo is not too close to 
the edges of the band, the contributions of the higher 
side bands become very small long before the energy Ej 
reaches the band edges. 

We will now calculate the different reflection and trans- 
mission amplitudes using Eqs. To first order in b n , 
only the first side bands with j = ±1 survive. If the 
oscillating potentials are given by 



^ S(x - x n ) b„ cos(wt + (f) n ) 



(8) 



we find that 



t R ,±i 



r R ,±i = 



I 

2vp~ 



t R {E ±1 ) b n< 



T'i'f 



x [ e Tml "/" F + rL (E )e~ i< ' k±1+k " +2kF ^ Xrl ] 



2v F 



t R (E ±1 ) t L (E ) 



Ti&n „ — i(k±i+k +2k F )x ri 



(9) 



Similar expressions can be derived for £l.±i and Tl,±i, 
but these will not be required below. The expressions in 
Eqs. © can be understood as arising from a sum over 
paths as explained in Ref. 37]. 

Given all the transmission and reflection amplitudes 
and 

right lead is given by [22|, [23j 



tR/L.j and r R /Lj, the^dc part of the current in, say, the 



N^l 2 {f R (E )-f R {E 3 )} 

+ \t Rij \ 2 {f L (E )-f R (Ej)}], 

(10) 

where q is the charge of the electron, and f a {E) = 



VI 



1] is the Fermi function in the lead a. At 



zero temperature, f a (E) = 1 if E < /j, a and if E > fi a . 
If we assume zero bias, [i R — fii = (we are defining en- 
ergy with respect to Ep), and restrict ourselves to terms 
of second order in the £>„'s, we obtain the expression 



L R,dc 



= q 



dEo 2 
dEo 2 

-7T- [\l-R,-l\ 
Z7T 



\^R,1 



\tR,-l\ 2 } 



(11) 



In the adiabatic limit uj — > 0, we can show from this 
that, up to order w, the current will involve only cross- 
terms of the form J2 m <n bmb n sin(0 m — 4> n ) multiplied 
by terms involving x m and x n . Let us call this the 'sin^' 
rule; it implies that two potentials cannot pump current 
if their phase difference is or n. We also see that if 
there is an oscillating potential at only site, there is no 
current at order uj. These features arise because of the 
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near cancellation in Eq. (TTT]) between \r R ,i 2 and 1 1"^, 1 1 2 

and between li^il 2 and |i_R,-i| 2 at order uj. However, the 
sm<j) rule does not hold, and an oscillating potential at 
even one site can pump current, if uj is larger than the 
resonance width; in that case Eq + uj and Eq — uj will 
be sufficiently different from each other so that (r^il 2 
and |rfl,_i| 2 or |t«,i| 2 and |t_R,_i| 2 are no longer almost 
equal to each other. The violation of the sin <j) rule will 
be explicitly demonstrated below for the case of pumping 
at two sites. 

We will now illustrate the above results for the case in 
which there are two static (^-function potentials given by 



a\5{x — x\) + a 2 5(x — x 2 ). 



(12) 



Defining it,; = cli/vf, we find the following expressions for 
the reflection and transmission amplitudes as functions 
of the wave number fc, 



= = 



1 



1 + i(ui + u 2 ) - u\u 2 {l - e i2 ( k+kp H^-zi)) ' 

i2(k+k F )x2 



r L = t R [-iu 1 j' z{h+kF)xi - iu 2 e 

+ u lU2 {e i2{k+kF)xi - e i2 ( fc +^)^)], 

r R = t R [~iu ie - l2( - k+kF)xi - iu 2 e- l2{k+kF)x2 

- u 1 u 2 (e- i ^ k+kF ^ xl - e -m+k F )x^ ( 13 ) 

These expression simplify considerably if u\ = u 2 = u, 
i.e., if we have a symmetric double barrier. Defining 9 = 
(k + kp)(x 2 — xi), we obtain 



\tR\ 2 = [ti| 5 



1 + Au 2 (cos 9 + usin( 



(14) 



We then see that \t R \ = |fj,| = 1, i.e., there is a resonance, 
whenever tan 6* = —1/u. This occurs at 9 = [n + 1/2)tt 
for u and at nir for u — > oo. For |u| 3> 1, the 
resonances are very sharp, with \t R \ 2 dropping to 1/2 
when one deviates from one of the resonant values of 9 
by l/(2u 2 ), namely, when kp deviates from one of the 
resonant values by l/(2u 2 \x 2 — xi\). 

Similar to other models considered in the earlier liter- 
ature, we can now study what happens when an electron 
is incident on the symmetric double barrier with an en- 
ergy Eq which is not at resonance, but the pumping fre- 
quency is such that Eq+uj or Eq — uj is equal to one of the 
resonant energies called E r . Then the oscillating poten- 
tials can change the energy of the electron to E r which 
can then transmit through the double barrier. Thus the 
pumping can help the electrons to access a resonance. 
We will compare below the cases of pumping at one site 
versus two sites, and also study the dependence of the 
pumped current on the phase difference <fi in the case of 
two-site pumping. 

In Fig. 1, we show the pumped dc current I Rt dc versus 
the Fermi energy Ep when there is a static and symmet- 
ric double barrier and an oscillating potential is applied 
at one site lying on the left of the double barrier. As Ep 
increases, we observe first a maximum at Ep = 0.883, 
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FIG. 1: (Color online) Current pumped through a symmetric 
double barrier by an oscillating potential applied at one site 
as a function of Ef- The barrier strengths m are both equal 
to 4, and their positions are at and ir. The resonant energy 
is E r = 0.917, and no bias is applied. The oscillating poten- 
tial is applied at the position — it with an amplitude 0.2 and 
frequency 0.035. We set vf = 1. The locations of the two 
maxima and the minimum in the middle are discussed in the 
text. 



then a minimum at 0.921, and then a small maximum 
at 0.960. These features arise for the followings reasons. 
The first maximum at 0.883 coincides with E r —ui within 
a spread given by the resonance width l/(27ra 2 ) ~ 0.010 
(note that this is significantly less than the pumping fre- 
quency uj — 0.035). This occurs because an electron ap- 
proaching from the left reservoir with an energy equal to 
Ep = E r — uj can get boosted up to the energy E r due 
to the oscillating potential; it can then transmit through 
the double barrier. The minimum at 0.921 coincides with 
E r ; note that this corresponds to a negative current, 
namely, the current is flowing to the left. This occurs 
due to a combination of two effects. An electron ap- 
proaching from the right reservoir with an energy equal 
to Ep = E r transmits through the double barrier; it can 
then get boosted to the energies E r ± uj by the oscillat- 
ing potential, which means that it cannot transmit back 
to the right through the double barrier. Similarly, an 
electron approaching from the left reservoir with the an 
energy equal to Ep — E r can get boosted to the energies 
E r ± uj by the oscillating potential, which means that 
it cannot transmit to the right through the double bar- 
rier. Hence, in both cases, the electron finds it easier to 
go to the left than to the right, leading to a net current 
to the left. Finally, the small maximum at 0.960 coin- 
cides with E r + uj. This maximum is not very robust; 
its height can change easily depending on the values of 
the various parameters because it is a result of several 
competing processes. An electron approaching from the 
right reservoir with an energy equal to E r (which is less 
than Ep) transmits through the double barrier; it can 
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then get boosted to the energies E r ±uj by the oscillating 
potential and then escape to the left reservoir. On the 
other hand, an electron approaching from the left reser- 
voir with an energy equal to E r ± u) can get boosted to 
the energy E r by the oscillating potential, which means 
that it can then transmit to the right through the double 
barrier. Depending on which of these is larger, the net 
current can be positive or negative. 




FIG. 2: (Color online) Current pumped through a symmetric 
double barrier by oscillating potentials applied at two sites. 
The barrier strengths are both equal to 4, and their positions 
are at and 7r. The resonant energy is E r — 0.917, and 
no bias is applied. The oscillating potentials are applied at 
the positions — tt and — 2n with the same amplitude 0.2 and 
frequency 0.035, but with a phase difference of n/2. We set 
Vf = 1. The locations of the maxima and the minimum are 
discussed in the text. 

In Fig. 2, we show the pumped dc current versus the 
Fermi energy Ep when there is a symmetric double bar- 
rier and oscillating potentials are applied at two sites 
(lying on the left of the double barrier) with a phase dif- 
ference of 7r/2 between the two potentials. Just as in 
Fig. 1, as Ep increases, we observe first a maximum 
at Ep = 0.895, then a minimum at 0.932, and then a 
small maximum at 0.963. The reasons for all these fea- 
tures are the same as the ones discussed above for Fig. 
1. We observe that the maximum value of the pumped 
current in Fig. 2 is about 10 times the corresponding 
value in Fig. I, showing that pumping by two oscillating 
potentials can be much more effective than by one oscil- 
lating potential. This can be qualitatively understood as 
follows. If oscillating potentials are applied to K sites, 
and their effects add up constructively, the reflection and 
transmission amplitudes in the first side band would be 
of the order of K times what one would get if there was 
pumping at only one site; this is evident from Eqs. ([§]). 
The pumped current would therefore be magnified by a 
factor of K 2 . 

In Fig. 3, we show the pumped dc current as a func- 
tion of the phase difference between oscillating potentials 
applied at two sites, for four different values of the Fermi 
energy Ep . These values have been chosen to lie between 
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FIG. 3: (Color online) Current pumped through a symmetric 
double barrier by oscillating potentials applied at two sites as 
a function of the phase difference (f>2 — (pi, for four different 
Fermi energies given by 0.893 (magenta solid line), 0.910 (blue 
dashed line), 0.920 (green dotted-dashed line) and 0.934 (red 
dotted line). The barrier strengths are both equal to 4, and 
their positions are at and tt. The resonant energy is E r = 
0.917, and no bias is applied. The oscillating potentials are 
applied at the positions — tt and — 2tt with the same amplitude 
0.2 and frequency 0.035. We set vf = 1. 



the first maximum and the minimum observed in Fig. 2. 
We observe that all the four curves are approximately of 
the form d\ + d,2 sin(</>2 — 4>\ + d^), where the values of 
di, g?2 an d ^3 are different for the different curves. All 
of these differ from a simple sin (f> rule which would cor- 
respond to the parameters d\ and d% being equal to 0. 
This violation of the smtp rule occurs here because we 
have chosen uj to be larger than the resonance width. 
The demonstration of this violation is one of the main 
results of this section. 

The results in this section can be generalized to finite 
temperatures by using the appropriate Fermi functions 
in Eq. (fl"U|) . One of the effects of a finite temperature 
T is to introduce a finite phase decoherence length given 
by l<p = hvp I (fcsT) [H]|. We expect this to destroy the 
resonance produced by the static double barrier when 1$ 
becomes smaller than the distance between the barriers; 
hence the peaks in the pumped current will become broad 
and eventually disappear. 



III. CHARGE PUMPING ON A FINITE RING 

In this section, we will study what happens when oscil- 
lating potentials are applied to a finite-sized ring which 
is not coupled to any reservoirs. This study may possibly 
be of interest in the context of transport in mesoscopic 
rings or even molecular rings fioj which are either not 
coupled to any reservoirs or are so weakly coupled to 
reservoirs that the momentum relaxation and phase de- 
coherence times are much longer than the time period of 
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the oscillating potentials. (A discussion of the possible 
experimental applications of our work will be presented 
in Sec. IV) . Once again, we will consider non-interacting 
spinless electrons. The Hamiltonian will be taken to be 
of the tight-binding form with nearest-neighbor hopping 
amplitudes and on-site potentials which are either static 
or oscillating with a time period T . Namely, we have 
H = H + V(t), where 

N N 

H = 7 (4 C «+1 + cl +1 C n ) + Y^ a n 4i°n, 

n—1 n—1 

N 

V(t) = ]T w n (t) c t c „, (15) 

n=l 

where w n (t) = w n (t+T) with T denoting the time period. 
(We will set the hopping amplitude 7 = 1 in all our 
numerical calculations). Note that we have included all 
the static potentials in Ho. We will impose the condition 
that 

/ dtw n (t) = (16) 
Jo 

for all n; if necessary, this can be ensured by adding 
a constant to w n (t) and subtracting the same constant 
from a n in Hq. We are therefore assuming that V satisfies 

f dt V(t) = 0. (17) 
Jo 

We will impose periodic boundary conditions, so that 
N + 1 = 1 in Eq. flTSj) . It will be convenient below to 
rewrite the operators in Eq. (fT5"j) as N x N matrices in 
the one-particle basis. Namely, 

(H )j k = - 7 (SjM+i +Sj t k-i) + a « 3hk$j,7i, 

TO 

(V(t)) jk = J2 hkhn- (is) 

TO 

We will now study various features of this model; in 
particular, we will calculate the current averaged over 
one time period for different choices of the oscillating po- 
tential V(t). We will assume that the system is not con- 
nected to any external reservoirs and has no mechanism 
for momentum relaxation or phase decoherence. We will 
begin with a given initial state in which the lowest p en- 
ergy levels of the Hamiltonian Hq in Eq. (fl"5"j) are filled 
(so that the filling fraction is p/N), and then evolve the 
system in time using only the total Hamiltonian H given 
by Eq. ([15)) . Consider the current operator for the bond 
(n,n + l), 

Jn = ~ i ( c n c n+l — C ri+l C ™)' (19) 

or, in matrix form, 

{Jn)jk = — i {§j,n$k,n+l — <5j>+l^fc,n)- (20) 



(Note that this matrix is Hermitian, imaginary and an- 
tisymmetric). We can then obtain the average value of 
the current at that bond by calculating the expectation 
value of J n over many time periods as described below; 
let us denote this average value by Idc- Note that in the 
absence of reservoirs, the system will not reach a steady 
state. The quantity l^ c should therefore not be thought 
of a steady state current; it is merely the current aver- 
aged over many time periods. We will see that apart 
from the average part called Idc, the current continues 
to oscillate in an aperiodic manner even if we wait for a 
period of time which is much longer than T. 

The time evolution of the system is governed by the 
unitary time evolution operator 

M 

U(t) = lim T TT e' lH ^ )d \ (21) 
i=i 

where T stands for time ordering (namely, 
TO(h)0(t 2 ) = 0(h)0(t 2 ) if tx > h), and we 
have divided the time t into M equal steps, i.e., 
dt = t/M and tj = (J — 1/2) dt; eventually, we have 
to take the limit M —¥ oo. We now consider U(T) 
where T is the time period. Let Vj and e l9j denote the 
eigenstates and corresponding eigenvalues of U(T), with 
j = 1, 2, • • • , N. Due to the periodicity of H in time, 
we see that U(sT)vj — e ls9j Vj for any positive integer s. 
Let us assume for simplicity that there is no degeneracy, 
so that e 4 ^ ^ e l9k if j ^ k; we can take the Vj's to form 
an orthonormal basis. 

We observe that if the periodic potentials in Eq. (fTS"]) 
are all shifted in time by the same amount, i.e., w n (t) —¥ 
w„(t + r), where r lies in the range [0,T], this gener- 
ally changes the unitary operator U(T). However, if we 
redefine the vectors Vj — > v'j — U(r)vj , we can use the pe- 
riodicity of w n (t) to show that the v'j will be eigenstates 
of the new operator U'(T) with the same eigenvalues e j . 
This can be proved as follows. Let us introduce a two- 
parameter notation for the time evolution operator 

M 

U 2 (t',t) = lim T TT e- iH ^ dt , (22) 

3=1 

where tj = t + (j - l/2)dt, and dt = (f - t)/M. Ac- 
cording to this notation, what we called U(t) earlier is 
actually U2(t,0) and the time shifted operator U'(T) 
is Ui^T + r, t). The periodicity of w n (t) then implies 
that U 2 {T+t 1 t)U 2 (t, 0) = U 2 (t, 0)U 2 (T, 0), namely, that 
U'{T)U(t) = U{t)U{T). Hence U{T)vj = e l0 i Vj implies 
that v'j = U{t)v satisfies U'{T)v' = e ld ^v' y Next, let us 
define the matrix for the current averaged over one time 
period, 

(J«)ik = \ J dtv)u\t)j n U{t)v k . (23) 

Then we can use the identity U(£)W{t) = U(t — r), for 
t > t, to show that 

v)u\t)J n U{t)v k = (v'^U f (t-T)j n U(t-r)v' k . (24) 
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Integrating both sides of Eq. (|24|) over one time period T, 
we see that the quantities (J n )jk defined in Eq. (|23|) are 
invariant under a shift in time by an arbitrary amount r. 

We now begin with an initial one-particle state ip a at 
time t = 0. This can be written as a linear superposition, 
ip a = Y^,j c ajVj , where c a j = Vjip a . If we evolve this initial 
state using the Hamiltonian H, the current averaged over 
an infinitely long time will be given by 

1 r sT 

JM = lim — / dtiplU\t)J n U{t)ip a 

s->oo si J Q 

= lim - V V e im ^- e ^c* aj c ak (J n ) jk . 

s— >oo s — — 

rn— 

(25) 

Since we have assumed that there is no degeneracy, i.e., 
lim - e m{ek ~ e ^ = 5 jtk . (26) 

S °° S m=0 

Hence Eq. (|23|) simplifies to 

N 

i=i 

We thus have a simple expression for the current averaged 
over an infinite amount of time, even though the current 
will oscillate in an aperiodic manner at all times due to 
the factor of e im ( e k-e 3 ) in Eq _ ^ Finally, if we start at 
t = with p electrons occupying the orthonormal states 
ip a , where a = 1, 2, • • • ,p, the current averaged over an 
infinite time will be given by 

p N 

J *= = E E m 2 ( J «)«- ( 28 ) 

o=l 3=1 

We have numerically computed the diagonal part of 
the averaged current {J n )jj for several different situa- 
tions. For the case in which H is real, we have made 
the following observations, all of which will be seen to be 
quite different from the situation on the infinite line, 
(i) If the periodic potentials w n (t) are shifted in time 
by an amount r as described above, the value of (J n )jj 
changes in general. This is because the overlaps \c a j\ 2 of 
the new eigenstates v'j with the initial states ipa is given 
by 

M 2 = \v)uHr)i, a \ 2 (29) 

which will generally depend on r, except in the trivial 
limit in which all the oscillating potentials w n {t) are set 
equal to zero. Hence the averaged current given in Eq. 
(|2"5| will vary with r, even though (J n )jj is independent 
of t as discussed earlier. All this is quite different from 
the situation on the infinite line where the time averaged 



current does not change if all the periodic potentials are 
shifted in time by the same amount. 

We note, however, that the dependence of Idc on the 
shift r (or, equivalently, the overall phase of all the os- 
cillating potentials) is weak if the amplitudes of all the 
oscillating potentials are small compared to the hopping 
amplitude 7. This statement will be made more precise 
below. 

(ii) If the periodic potential is time reversal invariant, i.e., 
if w n (t) — w n {~i) which is also equal to w n {T — t) due 
to the periodicity in time, then {J n )jj = for each value 
of j. We will present an analytical proof of this below. 
From this it follows that {J n )jj — even if the periodic 
potential is time reversal invariant only up to a shift, i.e., 
if w n (t) = w n {r — t) for some value of r. This too differs 
from the situation on the infinite line where the averaged 
current is non-zero even if the periodic potential is time 
reversal invariant, as long as the static potentials break 
the left-right symmetry. 

(iii) If the periodic potential is applied to only one site, 
then (J n )jj = even if the potential has no particular 
symmetry in time. We have found this numerically for a 
wide variety of potentials w n (t); (J n )jj always turns out 
to be zero to very high precision. We have no analyti- 
cal understanding of this remarkable result. This is also 
different from the situation on the infinite line where we 
have already seen that an oscillating potential applied to 
only one site can pump current if the static potentials 
break left-right symmetry. 

We will now prove analytically that (J n )jj = if H 
is real and w n (t) — w n (T — t) for all n. Following the 
notation of Eq. (I2T]) , this property of w n (t) implies that 

^ e -iH{t)dty = ( e -iH(T-t)dty^ ^ 

and therefore that U*(T) = U~ l {T). Hence U(T) Vj = 
e tBj Vj implies that U(T)v* = e l6j v*. The non-degeneracy 
of the eigenvalues of U(T) then implies that v* must be 
proportional to Vj. Hence we can multiply the Dj's by 
appropriate phases to ensure that v* = Vj for each value 
of j. Next, we can use Eq. (|30[) to show that 

U{T~t) = U*{t)U{T) (31) 

for any value of t lying in the range [0, T]. We can 
now combine Eq. (|31[) with the reality of Vj to show 
that vjW{T - t)j n U(T - t)vj, which must be real 
since J n is Hermitian, is equal to vjll^ (t)J*U(t)vj = 
—v'jU^(t)J n U(t)vj since J* = — J n . We thus see that 

(J n )j3 = i J dt v]u\t)J n U{t) V] = (32) 

due to an exact cancellation of the integrand at the times 
t and T - t. 

In Fig. 4, we show that the current integrated over one 
time period is zero if an oscillating potential is applied to 
only one site; as a specific example, we have considered 
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FIG. 4: (Color online) Instantaneous current at various times 
over one time period T — 10n (corresponding to ui — 0.2) for 
an oscillating potential applied at one site, for the third eigen- 
state of U(T) for a six-site system. Two cases are considered: 
the oscillating potential is (i) time-reversal invariant (blue 
dashed line), and (ii) not time-reversal invariant (red solid 
line). A static potential of strength 1 is applied at site 2, while 
the oscillating potential is applied at site 3. For case (i) with 
an oscillating potential cos(cjt) + 0.2cos(2u£), we see an exact 
cancellation between the values at times t and T — t, while 
case (ii) with an oscillating potential cos(oji) + 0.5sin(2o;t) 
does not show such a pair-wise cancellation. In both cases, 
the integrated current contributed by the third eigenstate is 
zero. 

the contribution to the current from the third eigenstate 
of U(T) for a ring with six sites, when a static potential 
and an oscillating potential are applied to two different 
sites. The way in which the current integrates to zero 
is quite different depending on whether the oscillating 
potential is time reversal invariant or not. If the oscil- 
lating potential is time reversal invariant, the integrated 
current vanishes due to a pair-wise cancellation from the 
times t and T — t as we have proved above. But if the os- 
cillating potential is not time reversal invariant, there is 
apparently no simple symmetry reason for the vanishing 
of the integrated current. 

Since the integrated current vanishes if an oscillating 
potential is applied to only one site, we will now study 
what happens when oscillating potentials are applied to 
two sites. In all the figures discussed below (Figs. 5- 
9), we will consider a six-site system in which a static 
potential is applied at site 2, and oscillating potentials 
of the forms bcos(u!t + <po) and bcos(ujt + tt/2 + <f>o) are 
applied at sites 3 and 4 respectively. 

We first study the dependence of the pumped current 
on the overall phase 0o of the oscillating potentials. We 
take b — 0.1, and use Eqs. (I28H29[) to compute the av- 
eraged current as a function of <p = 2ttt/T. In Fig. 5, 
we show the dependence of the current Idc on <pQ when 
the system has three electrons which corresponds to half- 
filling. We find that the variation of Idc from its mean 



value is about 8%. When we reduce b by a factor of 2, 
we find that the mean value of Idc decreases by a factor 
of 4 while its variation with <j)Q decreases by a factor of 
16. In the Appendix, we show that the mean value of 
Idc scales as b 2 which is in agreement with the numerical 
result quoted above. In addition, the numerics suggests 
that the variation of Idc with (j>o scales as b 4 . Thus, if the 
pumping amplitude b is small compared to the hopping 
amplitude 7, the variation of the pumped current with 
00 becomes much smaller than the current itself. 
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FIG. 5: (Color online) Total pumped current versus the over- 
all phase for oscillating potentials applied at two sites, for 
a six-site ring with three electrons. A static potential of 
strength 1 is applied at site 2, while the oscillating poten- 
tials are applied at sites 3 and 4 with an amplitude b = 0.1, 
a frequency uj = 0.2, and a phase difference of 7r/2. 

Next, we consider the contributions to the averaged 
pumped currents, Ij = {J n )jj of the different eigenstates 
of U(T), taking the overall phase to be given by 4>$ = 0. 
For the six-site ring, with a static potential of strength 
1 applied at site 2, the eigenvalues of the static part 
of the Hamiltonian, Ho, are equal to —1.8912, —1.0000, 
-0.7046, 1.0000, 1.3174, and 2.2784. In Fig. 6, we show 
the six contributions Ij as a function of the oscillation 
amplitude 6, when the oscillation frequency U) = 0.2 is 
non-resonant, i.e., it does not correspond to the difference 
between any two energy levels of Hq. We see that all the 
curves show a quadratic dependence on b, although two 
of them lie very close to zero; the quadratic dependence 
is in accordance with the results derived in the Appendix 
for the non-resonant case. In Fig. 7, we show the total 
pumped current Idc versus b when there are three elec- 
trons, where Idc has been calculated using Eq. (|28p and 
the individual pumped currents shown in Fig. 6. 

In Fig. 8, we show the contributions Ij as a func- 
tion of b, when the oscillation frequency cj = 1.7046 is 
resonant; it corresponds to the difference between the 
third and fourth energy levels of Hq. We see that four 
of the curves lie very close to zero and show a quadratic 
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FIG. 6: (Color online) Pumped currents for all the eigenstates 
of U(T) for oscillating potentials applied at two sites versus 
the oscillation amplitude b, for a six-site ring. A static po- 
tential of strength 1 is applied at site 2, while the oscillating 
potentials are applied at sites 3 and 4 with a phase difference 
of 7r/2 and a non-resonant frequency ui = 0.2. In all cases, 
the current is proportional to b 2 for small b. 
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FIG. 7: (Color online) Total pumped current for oscillating 
potentials applied at two sites versus the oscillation amplitude 
b, for a six-site ring with three electrons. All the parameters 
are the same as in Fig. 6. 



dependence on b, but the other two curves, which have 
large contributions from the third and fourth eigenstates 
of Hq, show a linear dependence on 6; the linear depen- 
dence agrees with the results derived in the Appendix for 
the resonant case. In Fig. 9, we show the total current 
Idc versus b when there are three electrons. This cur- 
rent varies quadratically with b for the following reasons. 
First, the two states which show a linear dependence on 
b in Fig. 8 are found to contribute with equal weight to 
Idc- Secondly, the slopes of these two curves in Fig. 8 



are equal and opposite at b = 0. Hence the linear depen- 
dences on b of these two curves cancel out when they are 
added up using Eq. (P251) to calculate Idc which there- 
fore shows a quadratic dependence on b. We also observe 
that the currents in Fig. 9 are about 90 times those in 
Fig. 7, for the same values of b; this shows the effective- 
ness of charge pumping at resonant versus non-resonant 
frequencies. 

0.02 | 1 1 1 1 1 1 1 1 i 1 
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FIG. 8: (Color online) Pumped currents for all the eigenstates 
of U(T) for oscillating potentials applied at two sites versus 
the oscillation amplitude b, for a six-site ring. A static po- 
tential of strength 1 is applied at site 2, while the oscillating 
potentials are applied at sites 3 and 4 with a phase difference 
of 7r/2 and a resonant frequency ui — 1.7046. In four cases, 
the current is proportional to b 2 , but in two cases, the current 
is proportional to b for small b (see text for details). 
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FIG. 9: (Color online) Total pumped current for oscillating 
potentials applied at two sites versus the oscillation amplitude 
b, for a six-site ring with three electrons. All the parameters 
are the same as in Fig. 8. 

Finally, let us briefly discuss what happens if the 
Hamiltonian is not real. For instance, if there is a mag- 
netic flux passing through the ring, its effects can be stud- 
ied by making the hopping amplitudes 7 complex in Eq. 
(fTijl) . In such cases, we find that there can be a net dc 
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current in the ground state even in the absence of any 
oscillating potentials; this is called a persistent current. 
The value of the dc current can then change even if we 
apply oscillating potentials which are time reversal in- 
variant or if an oscillating potential is applied at only 
one site. 



IV. DISCUSSION 

To summarize, we have studied charge pumping both 
on an infinite line and on a finite ring, for a system of 
non-interacting electrons at zero temperature. For the 
infinite line with reservoirs at the same chemical potential 
on the two sides, we have verified, in agreement with 
earlier work, that oscillating potentials applied to one or 
two sites can pump a net dc current as long as the left- 
right symmetry is broken by some static potentials; the 
oscillating potentials do not need to break time reversal 
invariance 22]. As a specific example, we have studied 
pumping in a situation where the Fermi energy differs 
from the resonant energy, but the pumping frequency is 
equal to the difference of the those two energies. We have 
contrasted the cases of pumping at one site and two sites. 
We have also shown that the sin (f> rule, which holds in the 
adiabatic limit (oj — > 0) for pumping by two oscillating 
potentials with a phase difference of <f>, fails if uj is larger 
than the resonance width. 

For a purely Hamiltonian evolution on a finite ring, we 
have shown that if both the static and time-dependent 
parts of the Hamiltonian are real, the net dc current is 
zero either if the oscillating potentials are time reversal 
invariant or if the oscillating potential is applied to only 
one site. We have presented an analytical proof of the 
former statement and have presented numerical evidence 
for the latter. It would be very useful if an analytical 
proof could be found for the vanishing of the current for 
one-site pumping with an arbitrary time dependence. 

We have argued that for a purely Hamiltonian evolu- 
tion, the system on a ring does not reach a steady state; 
the current averaged over one time period T continues 
to vary with time in an aperiodic manner although there 
is a simple expression for the current averaged over an 
infinitely long time. Further, the averaged current de- 
pends on the initial state and on the overall phase of the 
oscillating potentials since these two determine the over- 
lap between the initial state and the eigenstates of the 
evolution operator U(T). 

The facts that a steady state is not reached and that 
breaking of time reversal invariance of the oscillating po- 
tentials is required for charge pumping to occur on a 
finite ring but not on an infinite line is due to an im- 
portant difference between the models that we have as- 
sumed in the two cases. For the finite ring, we have 
assumed a purely Hamiltonian evolution, with no mech- 
anisms for momentum relaxation or phase decoherence. 
On the other hand, the study of the infinite line assumed 
that there are reservoirs which are maintained at certain 



chemical potentials. This implicitly assumes that there 
are mechanisms for energy or momentum relaxation; for 
instance, if an electron emerges from the double barrier 
with an energy equal to the resonant energy E r , we as- 
sume that when it reaches one of the reservoirs, it will 
relax down to the Fermi energy Ep if E r > Ep. It ap- 
pears that such relaxation processes, which we have not 
explicitly referred to in our calculations but which are 
necessarily present in the reservoirs, effectively lead to 
a breaking of time reversal invariance which is required 
to have a net dc current on the infinite line. An ear- 
lier study has shown that if a finite ring is coupled to 
reservoirs, then charge pumping can occur even if the 
oscillating potentials are time reversal invariant (2(| . 

It is also worth noting that on a finite ring, an electron 
which emerges to the right from the region containing the 
potentials (both static and oscillating) eventually comes 
back and enters the same region from the left. Thus the 
outgoing current on the right of that region must be equal 
to the incoming current on the left, when these currents 
are averaged over a long time. Similarly, the averaged 
outgoing current on the left of the region must be equal 
to the averaged incoming current on the right. These 
relations do not hold on the infinite line since an electron 
going out to the right (left) does not return to the left 
(right) of the region containing the potentials. 

We now turn to the possible experimental implications 
of our results. It is possible that the dependence of the 
averaged pumped current on the initial state and on the 
overall phase of the oscillating potentials on a ring with 
no mechanisms for momentum and phase relaxation can 
be observed experimentally. For mesoscopic rings which 
have a large number of impurities which lead to elas- 
tic scattering, the transport is diffusive and is character- 
ized by a diffusive round-trip time Td which is equal to 
3L 2 /(vpl c ), where l c is the mean free path and L is the 
circumference of the ring [65j . Another important length 
scale for ring systems is the phase coherence length 1$ 
which depends strongly on the temperature [6J, l66l . |67| . 
However, one may consider wires which have a sufficiently 
low density of scatterers and are at sufficiently low tem- 
peratures that both l c and 1$ are much larger L. In such 
a situation, we expect that if the pumping potentials are 
suddenly switched on, then the averaged pumped current 
will initially depend on the overall phase of the poten- 
tials. Eventually, at time scales which are much longer 
than l<p/vp and Td, the averaged pumped current will no 
longer remember the overall phase; however, a proof of 
this is beyond the scope of our analysis since we have not 
introduced any mechanisms for momentum and phase re- 
laxation. 

Another arena where our results could possibly be 
tested is the field of molecular electronics [42], l68| . We 
would like to propose an experimental set-up as follows. 
The conductance properties of aromatic molecules (which 
typically have a ring-like structure) have been studied 
extensively for several years. Typically, the transport 
properties of such molecules are studied by depositing 
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the molecule on a substrate and using an scanning tunnel 
microscope (STM) tip from above to probe the molecule 
[68[ . The number of electrons in the molecule can be fixed 
initially by bringing the STM tip close to the molecule 
and applying the correct potential to the tip. We pro- 
pose that the distance of the STM tip from the molecule 
can then be increased so that electron tunneling between 
the STM and the molecule becomes negligible; however 
the STM tip can still be used to induce an on-site pump- 
ing (oscillating) potential at a particular atom. This can 
be used to study the effect of pumping on the electronic 
transport through the molecule. A theoretical analysis 
of this would also require us to consider the effect of in- 
teractions between the electrons. 

An interesting problem for future studies may be to 
include the effects of interactions between the electrons 

69]. In particular, one can study whether interactions 
modify some of the peculiar features observed on the fi- 
nite ring, such as the absence of pumping for an oscillat- 
ing potential which is applied to only one site. In this 
context, we note that interactions between electrons are 
believed to play a role in determining the magnitude of 
the persistent current in rings placed in a magnetic field 

701. 
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Appendix 

We will use first order perturbation theory to derive ex- 
pressions for the eigenstates of U(T) on a finite ring and 
the integrated current in those eigenstates, for both the 
non-resonant and resonant cases. The treatment below 
will be seen to have interesting parallels with the use of 
first order perturbation theory to obtain the eigenstates 
of a time-independent Hamiltonian in the non-degenerate 
and degenerate cases respectively. 

The Hamiltonian of interest is H = Hq + V(t), where 
V(t) has a periodicity given by T = 2tt/uj. Let us first 
set V = 0. Hq will have a complete set of ortho normal 
eigenstates ipj and eigenvalues Ej, where we will assume, 
for simplicity, that the Ej's are non-degenerate. We will 
also assume that Hq is real; hence the ipj 's can be chosen 
to be real. If a unitary evolution operator Uq(T) is con- 
structed using Hq , its eigenstates and eigenvalues will be 
given by ipj and e~ lEjT . Note that since the wave func- 
tions ipj are real and J* = Jj = — J n , we have the useful 



relations 



ipjJni>j 

1p}j n 1pk 



iptJnipj for j ^ k. 



(33) 



We now turn on the time-dependent perturbation V 
which will be assumed to satisfy Eq. (fTTf . There are two 
different possibilities which we will discuss separately: 

(i) Non-resonant case where Ej — Ek is not an integer 
multiple of cj, i.e., e~ l ( Ej ~ Eh ' >T ^ 1, for any pair of states 
j, k, and 

(ii) Resonant case where Ej — Ek is an integer multiple 
of u), i.e., e~^ Ej ~ Ek ^ T = 1, but Ej / Ek, for some pair 
of states j, k. 

(iii) Resonant case where Ej = Ek for some pair of states 

Non-resonant case 

We will find the eigenstate of U(T), called Vi, which 
differs at first order in V from a particular eigenstate of 
Uo{T), say, ip\. We assume that Ui(t) = U{t)v\ has an 
expansion of the form 



N 
3=1 



iEit 



J ' 



(34) 



where we choose ci(0) = 1. (We are not worrying about 
the normalization of v\ here, although we see that v\ is 
normalized to 1 up to zero-th order in V). We expect 
that the deviation of c\ (t) from 1 at different values of t 
and also the values of Cj(t) for j ^ 1 will be of order V. 
From the Schrodinger equation 



i^p- = [H + V(t)} Wl (t), 

we find the following equations to first order in V, 
.d Cl (t) 



(35) 



dt 

dcAt) 
i , 
dt 



= ci(t) (iPtlVit)^), 
= d(t) (WTOlVi) e^- E ^ (36) 



for j ^ 1. At first order, we can replace C\{t) by ci(0) = 1 
in Eqs. (|36|) . This gives the solution 

ci(t) = 1 - i f 
Jo 

Cj (t) = aj - i [ dt'{iPj\V{t')\i) X ) e i ( E *- El)t ' ,{37) 
Jo 

for j ^ 1, where aj are constants of integration which 
can be fixed as follows. Since ci(T) = ci(0) due to Eq. 
([17)1. we see from the ip\ term in Eq. that «i(T) = 
e~ tElT Vi(0), i.e., v\ is an eigenstate of U(T) with eigen- 
value e~ lElT . We therefore demand that this should also 
be true for all the other terms ipj for j ^ 1 in Eq. (|34l) . 
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We therefore require that Cj(T)e~ 



iEjT 



e- tE ^ T c 3 (0) for 



all j ^ 1. This fixes the value of the constants ctj in Eq. 
(|37l) . and we find that 



Schrodinger equation in Eq. (|35[) now gives 
.dcx(t) 



Cj(t) = —i 



/ T dt'{^\V{t')\^) e *(^-Bi)f 



i(E 3 -Ei)T _ 1 



* / e^"^)*' (38) 



dt 

.dc 2 (t) 
1 dt 



d(t) (Vil^(i)IV'i) 
ca(t) (^|V(t)|^) 

+ d(t)e i05a - Bl >t (V 2 |V(t)^i). (42) 



for j ^ 1. Combining Eqs. ([51]). (|37|l and (|58"]>. we have 
an expression for the eigcnstatc of t/(T) to first order 
in V. Note that the corresponding eigenvalue remains 
g-iBiT £ Q -j-jjjg orc { er - this j s a consequence of the choice 
made in Eq. (fT7|) . 

We can now calculate the current averaged over one 
time period for one of these eigenstates, 



At first order, we can replace c\{t) and C2(t) by do and 
C20 on the right hand sides in Eqs. (|42|) . This gives 

ci(t) = do - i do I dt'i^lVit')^) 



u 20 / dt'^Vd')^) e^-^y 



(Jn)j 



1 

T 



ca(t) = c 20 - i c 20 / dt'{^ 2 ]V(i')l^2> 



dt v)u\t)J n U{t)v 3 . (39) 



ido / dt^WOhM ^-^''{43) 



Due to the fact that ^]j n ^ = 0, the current can only Uging Eq ^ we gee ^ 



get a contribution from cross-terms of the form tp^Jn^k 



for j ^ k arising from Eq. (|34j) . From Eqs. (|37H38I) . 
we see that (J n )jj has no contributions of order 1, while 
the contributions of order V can only come from a cross- 
term between ip\ and if>j for j ' 1. Such contributions 
are proportional to 



Cl (T) - d - ic 20 / dt<^i|V(t)|V>2> e^-^J*. 



C 2 (T) = c 20 -iq„ / <ft<Va|V(t)hM e« E *- E rt. 
Jo 

(44) 



(40) 



If we now demand that the state in Eq. (|41l) satisfies 
V (T) = e w v(0), i.e., that u(O) is an eigenstate of [7(T) 
with eigenvalue e l9 , we find that 



where Cj(t) is given in Eq. (l38|) . We can now do the 
integral in Eq. (j40|) explicitly and we find that it van- 
ishes. We therefore conclude that {J n )jj only receives 
contributions of second order and higher in V . 

Resonant case with E\ ^ E 2 

Let us now consider the case when two eigenstates 
of Uq(T), say, ipi and ip 2 , have the same eigenvalue 
e~ lT = e~ lE2T which implies that E% —E 2 is an integer 
multiple of w. We will assume, however, that E\ ^ E 2 . 
We will consider only the states 1 and 2, and will study 
how they can be combined to form eigenstates of U(T) 
to first order in V. Let us consider an expansion of the 
form 



d{T) e 



-iE±T 



c 2 {T) e 



-iE 2 T 



c 2 (0) 



d(0) 
This implies that 

- _ J T dt(^ 2 \V(t)\^) e *(«»-B0* 



C 2 Q 

do 



/„ dt(i/ji\V(t)\lp 2 ) ^(Bx-B^t 



Defining 



a = ™ 
T 



we see that 



dt{^\V{t)\^ 2 ) e 



i(E!-E 2 )t 



V(t) 



d(t) e' lElt Vi + c 2 (t) e-' ;B2 * V 2 , 



(41) 



C 2 Q 

do 



± 



(45) 



(46) 



(47) 



(48) 



where we now assume that both ci(t) and c 2 (t) are of 
order 1. (This is in contrast to the non- resonant case 
where only d(t) was taken to be order 1). Further, let 
us take d(0) — do an d c 2 (0) = c 2 o; we will assume 
that the deviations of d(t) an d d(t) from do an d do 
respectively will be of order V at all values of t. The 



so that c 20 /cio is a pure phase. Using Eqs. (|44M5|) . wc 
see that the ± sign in Eq. (|48|) corresponds to two solu- 
tions v±(t) which satisfy U(T)v±(0) = e l0± v±(O), where 



e i0 ± = e -,E lT 



i ■ rr °20 

1 — lOtl — 

do 



-iT (E 1 ±\a\) 



(49) 
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to first order in a. We thus see that the degeneracy of 
eigenvalues of U (T) is broken at first order in V, with 
the phases 9± being split by equal and opposite amounts. 

Finally, we can compute the current averaged over one 
time period as defined in Eq. (|39[) for either one of the 
states, v+ or u_. Once again, only cross-terms of the 
form tp\ j n ip2 will contribute. We find that terms of order 
1 vanish because they are of the form 



(50) 

and J Q fie'' £l_fi2 '' = because E\ — E 2 is an integer 
multiple of uj but E\ ^ E 2 . However, unlike the non- 
degenerate case, there is now no reason for contributions 
of order V to vanish in general. Hence (J n )jj can get 
contributions at first order in V. 



ci(T) = ci(0) = Cio and c 2 (T) = c 2 (0) = c 20 due to Eq. 
(|17|) . Further, there is now no relation between cio and 
C20; we can choose cio and c 2 o in an arbitrary way to ob- 
tain two orthonormal states v± which are eigenstates of 
U(T) with the same eigenvalue e~ lElT . Thus the eigen- 
values of U(T) and their degeneracy do not change to 
first order in V. 

When we compute the current averaged over one time 
period as denned in Eq. (|39|) for either one of the states, 
i> + or we find that, depending on the choices of c 10 
and c 2 o, we can get a contribution to zero-th order in V. 
To be explicit, this is given by 



dt [c* c 20 ip{j n ip2 + c 20 ci V4*Wi 
T (c* 10 c 20 - c 20 ci ) tp\j n ^2, 



(51) 



Resonant case with E\ = E 2 

Finally, let us consider the case when two eigenstates 
of Hq, say, ipi and ip 2 , have E\ = E 2 . This implies that 
they also have the same eigenvalue e~ tElT — e~ lE2T of 
Uq(T). As in the previous subsection, we will study how 
these two states can be combined to form eigenstates of 
U(T) to first order in V. 

We can check that all the discussion from Eq. (T4"Tj) to 
(1441) will remain valid in this case, except that we have 
to substitute E\ — E 2 everywhere. We then find that 



where we have used Eq. (|33|) . 

It may seem strange that if E\ — E 2 , there can be 
a non-zero current in the states v + and t>_ separately 
even at the zero-th order in V. However, we note that 
an energy degeneracy usually does not occur if static po- 
tentials are present. But if there are no static potentials 
present, then the system is translation invariant, and we 
see that there are momentum eigenstates (with momenta 
±fc) which carry equal and opposite currents even in the 
limit that the amplitudes of the oscillating potentials go 
to zero. 
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